Reliable Electronic and Magnetic Structure for BaFe 2 As 2 using the PAW method 
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By carefully choosing parameters and including more semi-core orbitals as valence electrons, we 
constructed a high quality projected augmented wave (PAW) dataset that yields results comparable 
to existing full-potential linearized augmented plane- wave calculations. The dataset was then applied 
to BaFe2 As2 to study the effect of different levels of structure optimization, as well as different choices 
of exchange-correlation functionals. It is found that the LDA exchange-correlation functional fails 
to find the correct SDW-AFM ground state under full optimization, while PBE obtains the correct 
state but significantly overestimates the magnetism. The electronic structure of the SDW-AFM 
state is not very sensitive to structure optimizations with the PBE exchange-correlation functional 
because the position of the As atoms are preserved under optimizations. 



The recent discovery of iron pnictide superconductor 
parental materials LnFeAsO (1111)±, ^eFe 2 As 2 (122)2, 
and ^4FeAs (lll)&i has stimulated tremendous inter- 
est and attention starting about a year ago. From the 
very beginning, various efforts have been made in order 
to determine the structure and electronic and magnetic 
properties of these materials from first principles". 
Both the local density approximation (LDA) and gener- 
alized gradient approximation (GGA) to the exchange- 
correlation functional have been employed, while for 
electron-iron interactions, various methods including 
ultrasoft pseudopotential£&2 (USPP), projected aug- 
mented waveiS (PAW), and full-potential linearized aug- 
mented plane- wav e 5 ' 7 ' 11 (FLAPW) were used. While 
FLAPW is generally regarded as the most accurate 
method among the three, it is computationally much 
more expensive than both USPP and PAW methods. 
Therefore, study using FLAPW often does not involve 
full relaxation including the unit cell parameters. For 
USPP and PAW, previous studies showed that both 
methods tend to overestimate the local magnetic moment 
on individual Fe atoms as well as the relative energy dif- 
ferences between different magnetic states^', indicating 
a serious problem regarding the commonly-used current 
pseudopotential or PAW datasets. Although it is now 
widely perceived that the calculated electronic structure 
is only comparable with experiments when the crystal 
structure is fixed as experimentally measured, this state- 
ment has never been verified under full relaxation using 
FLAPW or reliable USPP/PAW calculations. It is also a 
question for this system if the LDA exchange-correlation 
functional is superior to GGA or vice versa. In this pa- 
per, we report our efforts to generate a reliable PAW 
dataset for these materials, exemplified by BaFe2As2- 
We compare all of our results with FLAPW results in 
detail whenever feasible, and we study the influence due 
to full relaxation of the crystal structure as well as dif- 
ferent exchange-correlation functionals on the calculated 
results. 

In principle, one can approach all-electron accuracy 
with pseudopotential or PAW methods by including more 
semi-core orbitals into valence orbitals, but in practice it 
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USPP 


atom 


valence 
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Fe 
Ba 
As 


3s 3p 3d 4s 4p 
5s 5p 5d 
4s 4p 3d 


2.0 (2.1 d ) 
2.9 
2.1 


3d 4s 4p 
5s 5p 6s 
4s 4p 


2.2 (2.3 d ) 
2.5 (2.P) 
2.4 (2.1 p ) 



TABLE I: parameters used to generate the PAW datasets in 
this paper. r cut unit is in bohr, the number in the bracket 
indicates the value used for the orbitals with angular momen- 
tum represented by the superscript. 



is prohibitively difficult to do so for the ultrasoft pseu- 
dopotential method. Nevertheless, this process is prac- 
tical with Bloch's PAW method^. For this study, we 
follow the slightly modified recipe by Kresse et al*£, as 
implemented in the PWSCF package^. The valence or- 
bitals (including the semi-core orbitals) included in the 
datasets are listed in TABLE U compared with previous 
pseudopotentials. The generated dataset was carefully 
tested with logarithmic derivatives, ionization potentials, 
and simple crystal structures. The logarithmic derivative 
energy range considered was taken to be from —7.0 Ry to 
3.0 Ry, and orbital angular momentum I from to 4 (s, 
p, d, f, g). These tests ensured the quality of the PAW 
dataset. Throughout this paper, we used a cutoff en- 
ergy of 40 Ry for plane- waves in both PAW and USPP 
calculations to ensure convergence, and the irreducible 
Brillouin zones were sampled with the Monkhorst-Pack 
scheme^, 8x8x8 for the body-centered tetragonal 
(bet) unit cell (non-magnetic NM/checker-board anti- 
ferromagnetic CB-AFM states), 6x6x8 for the base- 
centered orthorhombic (bco) unit cell (spin-density wave 
antiferromagnetic SDW-AFM state). Both LDA and 
Perdew-Burke-Ernzerhof (PBE) parameterization*^ to 
GGA were used. 

In order to directly compare with existing literatures 
using FLAPW method, we calculated the BaFe2As2 
parental compound with different structural parameters. 
These calculations were performed on the structures used 
in Ref. [lg, using NM lattice parameters with relaxed za b 
and with experimental z\ s . It should be noted that both 
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FLAPW 


1.75 


0.70 


1.60 


N/A 


PAW 


1.75(2.34) 


0.60(1.60) 


1.54 


N/A 


USPP 


2.12 


1.04 


1.90 


0.58 


AE (meV) 










FLAPW 


92 


N/R 


41 


N/A 


PAW 


65(174) 


< 1(28) 


19 


N/A 


USPP 


119 


6.9 


60 


< 1 



TABLE II: Comparison between our PAW-LDA results with 
FLAPW 10 and USPP results with the same structure used in 
Ref. [H. zas refers to optimized zas with experimental lattice 
constants, and full-opt means fully optimized structure with 
relaxed lattice constants. AE is the total energy relative to 
the NM state per iron atom. The PBE results are given in 
parentheses. 
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FIG. 1: Comparison between FLAPW band structure and our 
PAW result for the non-magnetic state Bal22 with relaxed za s 
but fixed experimental lattice parameters (LDA result). The 
PAW result (red solid line) is almost exactly the same as the 
FLAPW result (blue dashed line). 

these parameters were used for the SDW-AFM state as 
well. Results are listed in Table [TTJ Our PAW calcu- 
lations agree with FLAPW results very well including 
the local magnetic moment on Fe atoms, whereas USPP 
always exagerate the magnetic moment, as reported pre- 
viously. Using the LDA version of these datasets, we 
also calculated the band structure for the NM configu- 
ration relaxed- zas structure. We compare our result and 
the FLAPW result in Fig. [H and it is apparent from 
this plot that our datasets almost exactly reproduced 
the FLAPW band structure. From Table HH it is also 
observed that using the correct experimental lattice pa- 
rameters for SDW-AFM configuration will remove some 
stress-induced magnetization on Fe atoms, and a further 
relaxation of za s "will eventually bring the local magnetic 
moment down to 0.4/Lts/Fe. In this situation, the SDW- 
AFM state energy is lower than the NM state energy by 





NM 


SDW-AFM 




Opt-ZAs full-Opt 


Opt-ZAs 


full- opt 


a 


3.963 3.873(3.968) 


5.615 


(5.697) 


b 


3.963 3.873(3.968) 


5.574 


(5.594) 


c 


13.02 12.14(12.49) 


12.94 


(12.76) 


ZAs 


0.341(0.344) 0.347(0.346) 


0.341(0.352) 


(0.3514) 


m Fc (^s) 


0.0 0.0 


1.62(2.14) 


(2.19) 


AE (meV) 


0.0 0.0 


105 (83) 


(68) 



TABLE III: Physical properties after different levels of struc- 
tural optimization using PAW. The SDW-AFM state does not 
exist after full relaxation using LDA, thus only PBE results 
are provided. The values in parentheses are obtained with 
PBE. 



22 meV/Fe. Neither FLAPW nor PAW is able to obtain 
a CB-AFM configuration with the relaxed za s , but this 
configuration remains stable using USPP. 

Next we evaluate different exchange-correlation func- 
tionals and study the band structure dependency on the 
degree structural optimization. As summarized in Ta- 
ble IIII1 the PBE exchange-correlation functional always 
exagerates the local magnetic moments on Fe atoms, 
and the magnetic moment remains relatively unchanged 
with respect to the structural optimizations. However, 
it should be noted that the PBE exchange-correlation 
functional yields crystal structures much closer to exper- 
iment, and the za s parameters are significantly differ- 
ent using PBE and using LDA. PBE also yields a large 
magnetic moment change upon zas optimization when 
the LDA crystal structure is used; see Table |TT] data in 
parentheses. It was also observed that a full optimization 
of both lattice parameters and za s within LDA fails to 
yield a SDW-AFM ground state (in fact, it does not show 
any magnetic instability after full optimization) due to 
the fact that LDA underestimates bond lengths/lattice 
constants while this system is too delicate to allow such 
underestimation. PBE overestimates both parameters, 
obtaining the correct ground state of the system under 
full optimization due to an error cancellation. Never- 
theless, both exchange-correlation functionals underesti- 
mate the c-axis in both the NM and SDW-AFM states, 
indicating an incompleteness of LDA and GGA exchange- 
correlation functionals. Since the local magnetic mo- 
ment of Fe atoms depends strongly on za s parameters, as 
shown in previous studies and in this work, we conclude 
that for the magnetic properties, neither LDA nor PBE 
is capable of describing the subtleness of this particular 
system; but PBE correctly describes the ground state of 
the system under full optimization. 

Due to the loss of the SDW-AFM ground state in LDA 
after full optimization, the following discussion includes 
PBE results only. As pointed out by Ref. [lj]and[ll, DFT 
calculations can be comparable with experimental results 
only if zas is fixed at the experimental value. Thus, these 
calculations were performed with FLAPW without full 
optimization of lattice parameters. Figure [2] depicts a 
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FIG. 2: Band structure for NM state Bal22 with (a) exper- 
imental structure, and (b) fully optimized lattice parameters 
as well as za s using PBE exchange-correlation functional. The 
band structure for the opt-ZAs structure looks very close to the 
LDA result (Fig. [I]), and is therefore not shown here. The 
effect of structural optimization mostly appears in the thicker 
blue band. 



set of band structures corresponding to different levels of 
structural optimization. Apparently, optimization does 
not alter significantly the electronic structure around the 
M point, where two electron pockets were found. More- 
over, even the size of these electron pockets does not 
strongly depend on the structural optimization. In con- 
trast, the structural relaxation substantially changed the 
band structure around T and Z. The main effect of struc- 
tural relaxation is optimization pushes one of the Fe 3d- 
bands toward lower energy (the thicker blue solid line in 
Fig. [2]); thus, only the two hole pockets from the rest of 
the 3d-bands remains after za s optimization as well as 
after full structural relaxation. 

Unlike the NM state, the band structure of SDW-AFM 
state does not exhibit large variation before and after 
structural optimization (Fig. [3]). This can be easily un- 
derstood since za s in SDW-AFM state was maintained 
during structure relaxation. We have used the SDW- 
AFM experimental structure and calculated the band 
structure with LDA, and there is not much difference 



FIG. 3: Band structure for SDW-AFM state Bal22 with (a) 
experimental structure, and (b) fully optimized lattice pa- 
rameters as well as za s using PBE exchange-correlation func- 
tional. Due to the spin degeneracy, only the a spin band 
structure is shown here. 



between the LDA result and the ones shown in Fig [3l 
despite of the significant difference (0.5 /^s) in the tofo 
found in these calculations. It is then apparent that the 
fundamental discrepancy between the LDA and the GGA 
calculations in the previous literature is due to the dif- 
ference between the structure optimized with LDA and 
GGA, respectively. The full optimization with the PBE 
exchange-correlation functional does not significantly al- 
ter the band structure in the SDW-AFM state, but is 
crucial to study realistic doping and pressure effects in 
order to be consistent. 

Finally, we present the Fermi surfaces reconstructed 
using maximally localized Wannier function o 19 ' 20 (ML- 
WFs) (Fig. 2]). For the NM state, as we have seen in Fig. 
[2l the optimization leads to very different band struc- 
tures, thus results in very different Fermi surfaces as well. 
After full-optimization, the size of hole pockets around T 
is greatly reduced and the two electron pockets are much 
less affected, as suggested by the band structure calcula- 
tion, and therefore the expected (7r,7r) Fermi-surface nest- 
ing is missing. For the SDW state, since the optimization 
does not significantly alter the electronic structure as it 
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FIG. 4: Fermi surfaces reconstructed using maximally local- 
ized Wannier functions (MLWFs). k x and ky denote the re- 
ciprocal vectors corresponding to a and b for the conventional 
tetragonal unit cell. 



does in the NM state, we thus show only the Fermi sur- 
faces after full optimization. Unlike its NM state, the 
SDW Fermi surfaces are surprisingly simple and do not 
show very interesting features such as Fermi-surface nest- 
ing. Two extended electron pockets were found, and two 
hole pockets appear around k z = 0.5. 



In conclusion, we have constructed a reliable PAW 
dataset that yields results comparable to the FLAPW 
method, that is then used to study electronic and mag- 
netic structure for BaFe2As2- A thorough comparison 
between exchange-correlation functional as well as struc- 
tural optimizations were presented. It is found that the 
LDA exchange-correlation functional will miss the SDW- 
AFM configuration under full structural optimization, 
and the PBE gives the correct ground state. The im- 
portant c-axis lattice parameter and zas can be obtained 
within reasonable error for the SDW-AFM state with the 
PBE exchange-correlation functional, but both values are 
noticeably reduced for the NM state. Therefore, the elec- 
tronic structure of the SDW-AFM state is insensitive to 
structural optimization, but the optimization will yield 
unacceptable electronic structure for the NM state. PBE 
can describe the SDW-AFM ground state very well but 
not the NM state, and the LDA fails for both. Based on 
the apparent correlation between magnetism and crys- 
tal structure, we therefore propose that LDA/GGA fails 
for the NM state because of a strong magnetic-phonon 
coupling within the iron-pnictidc compounds, and that 
this strong coupling may be the mechanism for the elec- 
tron pairing. We suggest future study to focus on the 
magnetic-phonon coupling effect. 
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